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l«  ah-,  1  m\C  t 

‘A  new  theory  of  turbulent  shear  flow  is  based  on  the  discovery  of  a  second  boundary 
layer,  which  we  call  the  viscous  mesolayer*  of  thickness  proportional  to  Taylor's  micro¬ 
scale  for  tl^  outer  flow  or,  scaling  on  friction  velocity  ut  and  viscosity  v,  proportional  to 
Ra  =  (u  a/v)‘,  where  a  is  the  outer  length.  *rhc  existence  of  this  layer  between  the  core  and 
the  subfayer  makes  it  impossible  to  match  the  outer  (core)  expressions  for  mean  quantities 
with  the  inner  expressions  based  on  Prandtl's  "law  of  the  wall",  as  in  classical  theory. 
Indeed,  the  supposed  region  of  applicability  of  the  classical  theory,  namely  the  layer  in 
which  both  inner  and  outer  behaviors  are  imposed,  is  precisely  in  the  mesolayer  where 
neither  behavior  may  be  imposed. 

The  mesolayer  concept  is  used  to  develop  a  new  theory  of  the  distribution  of  mean 
quantities  in  pipe  flow.  Solutions  are  obtained  for  the  mean  quantities  in  an  inertial  sub¬ 
range  and  in  a  viscous  subrange  on  cither  side  of  the  mesolayer.  Major  alterations  in  the 
classical  theory  result.  x 
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A  NFW  TIIFORY  OF  ITRRl’l.FNT  Fl-OW  IN  A  PIPF 


AUSTllUT 


A  now  t heory  of turbulont  shear  flow  is  based  on  (ho  ilisoovory  of  a  second 

boundary  la.vor,  whioh  wo  oall  tho  viscous  mosolayor,  of  fhioknoss  proportional  (o 

Taylor's  miorosoalo  for  tho  outer  flow  or,  sealing  on  friotion  volooity  u  ami 

\  ,  1 

viscosity  v,  proportional  to  R“  (U  a  vP,  whoro  a  is  tho  out  or  length.  Tho  oxistonoo 
ot  this  lavor  botwoon  tho  ooro  anti  tho  sublayor  makes  it  impossiblo  to  matoh  tho 
°*it or  (oorol  oxprossions  for  moan  tpiantitios  with  tho  i tutor  oxprossions  hasotl  on 
Prandtls  "law  of  tho  wall",  as  in  olassioal  theory.  Indootl,  tho  supposed  rogion 
of  applicability  ot  tho  olassioal  thoory,  namoly  tho  lavor  in  whioh  both  innor  anil 
out  or  bohaviors  aro  imposotl,  is  precisely  in  tho  mosolayor  whoro  noithor  bolwvivu’ 
may  bo  impasotl. 

Faoh  moan  quantity  is  oxpaitilori  in  two  infinito  sorios,  tho  first  forms  of 
whioh  aro  appropriato  for  tho  out  or  rogion  whioh  is  now  .•  ss  It*  anti  for  tho  wall 
rogion,  0(11,  whoro  :  is  a  now  similarity  variablo,  7  R*  ,  anil  .•  is  soaloil  on 
u  and  v.  Fquatingtho  two  sorios  solvos  for  oaoh  of  tho  moan  quant  it  ios  subjoot 
to  ohoioos  of  univorsal  oonstants.  It  is  possiblo  to  ohooso  thoso  oonstants  to 
ponnit  oloso  agroomont  with  tho  data  ovor  most  of  tho  flow.  Additional  ovidonoo 
for  tho  oxistonoo  and  importation  of  tho  mosolayor  may  b<>  found  in  tho  data  for 
pipo  flow  in  whioh  (11  thoro  aro  two  ix'aks  of  turbulotioo  onorg>-  in  tho  wall  rogion, 
(oorrosponding  to  a  first  and  sooond  oomponont  of  turbulont  motionl  as 
prodiotod  by  thoory;  (.’1  tho  rms  normal  volooity  tT«  (dominatod  in  tho  thoory  by  tho 
first  oomponont  ovor  a  fhiok  rogion  noar  tho  walll  is  Roynolds  nnmbor  indoix'ixlont  in 
a  lavor  thick  compared  to  tho  regions  in  whioh  tho  first  oomponont  of  a,,  and 


i. 


ov  dominates,  as  predicted  by  theory;  (111  oN  and  o,  hot h  reach  maxima  at  or  near 

;  v 

z  a  H,  ,  z  a,R*f  ,  and  a?  is  considerably  less  than  a, ,  all  of  which  are  predicted 
by  theory;  (41  o„  (dominated  by  the  first  component  in  a  thin  layer  near  the  walll 
rises  to  a  maximum  anti  then  decreases,  seemingly  approaching  a  constant.  It  then 
begins  to  increase  again  (because  of  the  second  component!.  The  curve  lias  an  in¬ 
flection  |H'int  and  then  another  maximum  ami  this  inflection  i>oint  occurs 
at  z  a* R*  as  predicted;  t ;»>  tht'  Reynolds  stress  peaks  at  z  a4  R“  as  predicted; 

(M  tht>  theory  predicts  that  there  is  a  region  just  above  the  sublayer  in  which  Ti 
increases  as  l’  ,  in  z.  At  a  (xont  in  tht'  mesolayer  there  is  a  change  to  a  smaller 
slope,  (t-  ,1  H  ,h'1  <n  z  where  is  negative,  ami  this  is  followed  by  increasing 
slopes  in  the  outer  region.  This  complicated  pattern  is  a  weak  variation  about 
the  classical  in  i  behavior  and,  of  course,  is  unattainable  in  classical  theory, 
but  the  predicted  behavior  appears  in  the  data. 

Many  rentiers  will  feel  reluctant  to  abandon  the  classical  theory  even  in  the 
face  of  the  success  of  the  new  theory  in  predicting  the  behavior  of  many  measured 

mean  quantities.  The  paper  contains  a  discussion  indicating  that  the  location  of 

' 

the  peak  of  the  Reynolds  stress  at  z  ■»-  R‘  directly  contradicts  one  of  the  basic 
assumptions  of  classical  theory  regarding  the  nature  of  the  "buffer  layer"  between 
inner  and  outer  layers. 

The  paper  contains  a  conjecture  that  the  "bursts"  are  characteristic  of  the 
mesolayer  and  tliat  they  grow  as  they  move  outward  and  along  tht'  wall.  Indeed, 
if  they  grow  by  molecular  diffusion,  they  grow  to  a  size  of  the  order  of  the  meso¬ 
layer  thickness  in  the  time  period  of  a  large  core  eddy.  Finally,  an  argument  is 

_  i. 

advanced  that  the  intermitteney  >  is  proportional  to  R  ■*,  yielding  v  -  0.  I  -  0.  J, 
at  typical  R  ,  in  rough  agreement  with  observations. 
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A  NEW  THEORY  OF  TURBULENT  FLOW  IN  A  PIPE 

1.  The  Problem  of  Turbulent  Shear  Flow. 

The  classical  theory  of  turbulent  shear  flow  near  a  surface  (Ka'rman,  1930, 
Prandtl,  1932,  Monin  and  Yaglom,  1971,  Chp.  3)  is  enduring  (despite  the  clear 
experimental  evidence  that  it  is  too  simplistic'  because  nothing  more  appealing 
has  been  offered  and  because  the  theoretical  mean  velocity  and  the  drag  co¬ 
efficient  in  a  pipe  seem  close  to  observations.  Its  predictions  for  the  root- 
mean-square  velocities  are  not  as  good  and  recent  efforts  have  been  made  by 
Townsend  (1976)  and  Perry  and  Abell  (1977)  to  improve  the  theory  for  these. 
Many  other  efforts  have  been  made  to  develop  higher-order  approximations  for 
various  mean  quantities  assuming  that  the  classical  theory  is  a  correct  first 
approximation,  for  example  Afzal  and  Yajnik  (1973).  Most  such  investigations, 
including  those  of  the  present  paper,  lead  to  mathematical  expressions  involving 
one  or  more  universal,  unknown  constants,  and  it  is  usually  not  convincing  to 
demonstrate  agreement  with  the  data  over  the  limited  ranges  of  the  theories 
when  the  constants  are,  in  fact,  chosen  to  force  agreement  at  one  or  more 
points  ot  the  data  curves.  The  present  paper  contains  a  theory  which  hopefully 
passes  more  severe  tests. 

A  useful  approach  to  the  classical  theory  was  developed  independently 
by  Izakson  (1937)  and  Millikan  (1938)  who  assumed  two  regions  of  flow,  an 
"inner"  region  near  the  surface,  and  an  "cuter"  region  further  out.  In  the 
inner  region  it  is  assumed  that  the  outer  length  parameter  a  (e.  g. ,  the 
radius  of  a  pipe)  is  unimportant  for  all  mean  quantities.  This  assump- 
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tion  was  first  advanced  by  Prandtl  (1925,  1932)  and  is  called  the  "universal  law 
of  the  wall”.  We  call  this  assumption  "independence  of  outer  length".  In  the 
outer  region,  it  is  assumed  that  the  viscosity  coefficient  v  is  unimportant  for  mean 
velocity  shears  and  for  the  velocity  defect  (Uo-u)  and,  indeed,  for  all  mean 
quantities  not  connected  directly  to  the  small  scales  of  the  turbulence.  This  is 
called  "Reynolds  number  similarity"  (Townsend,  1976)  anil  was  also  formulated 
first  by  von  Karman  (1930).  We  call  this  assumption  "independence  of  vis¬ 
cosity".  Roth  assumptions  have  been  verified  bv  measurements,  for  example 
l-  ritsch  (192S),  Laufer  (1954),  and  Perry  and  Abell  (1975).  A  third  assumption 
is  that  there  is  a  finite  region  of  overlap  of  the  two  regions.  We  discuss  this 
assumption  at  length  in  Section  4  and  find  it  wanting.  Nevertheless, 
using  these  three  assumptions,  mathematical  operations  similar  to  those  used 
in  Section  3  lead  to  the  fundamental  results  for  the  overlap  region,  namely 


u 

-  fn  /.  +  A 

(i) 

H 

- 1) 

--  in£  +  R 

X  v 

(2) 

Uo 

-  in  R  +  A+  R 

X  T 

(3) 

1.  o„ 

Al*  °v  •V'»  lTw  A3 

(4) 

where  <ru  ,  irv ,  ir„  are  root -mean- square  velocities,  fis  distance  from  the  wall, 
scaled  on  a,  and  all  other  quantities  are  scaled  on  the  inner  variables  v  and  friction 
velocity  at  the  wall  u  .  The  Reynolds  number  used  here  is  uTa/v  Rt  and  u0  is 
the  non-dimensional  outer  velocity  which,  for  pipe  flow,  is  the  velocity  at  the 


L 


center  lino,  l'he  constant  .is  called  von  barman's  constant  am!,  except  for  one  set 
of  measurements  (Basinger,  et  al. .  l!)7l)  .  appears  to  be  about  the  same.  '  O.  -U. 
in  al!  laboratory  and  geophysical  flows. 

The  agreement  of  (1  >-(3)  with  the  data  is  puzzling  because  close  observation 
of  phenomena  in  shear  flow  over  the  past  15-20  years  indicates  tint  the  structure  of 
the  turbulence  is  very  different  from  the  picture  inherent  in  the  classical  theory. 

T°  bo  specific,  all  the  results  in  (l)-(-H  may  be  obtained  bv  dimensional  analysis 
assuming  independence  ot  both  .  and  a  in  l  •  ■  z  ■  -  it  and  this  approach  predicts 
turther  t  hat  all  length  scales  in  the  overlap  region  are  of  the  on  let-  of  the  distance  z 
from  the  surface.  Measurements,  however.  (Mitchell  and  llanratty,  Ukm,  Tritton, 
Ith'.i  indicate  that  the  longitudinal  length  scale  is  much  larger  than  z  near  the 
surface  and  may  be  of  the  order  of  the  outer  length  scale  a.  We  know  for  certain 
now  that  large  eddies,  associated  with  motions  in  the  core,  are  exerting  a  strong 
influence  in  the  region  just  above  and  even  in  the  viscous  sublaver.  Another 
example  is  the  appearance  of  the  "bursts"  above  the  viscous  sublayer  which  have 
a  time  period  much  longer  than  the  only  time  period  of  the  classical  theory,  namely 
distance  from  the  surface  divided  by  friction  velocity,  and  which  seems  to  scale  on 
the  outer  variables  (Laufcr  and  Narayanan,  i‘)7i,  too,  et  al. .  Obviously, 

then  the  independence  of  a  assumed  in  the  so-called  buffer  region  is  untenable. 

The  existence  and  importance  of  large  eddies  near  a  wall  have  boon  demon 
strated  by  theoretical  and  experimental  studies  of  the  distortion  of  turbulence  in 
a  wind  funnel  by  a  surface  parallel  to  the  flow  and  moving  at  the  speed  of  the  free 
stream  to  yield  shear  free  turbulence  (Uzkan  and  Reynolds,  i:h>7,  Thomas  and 
Hancock,  197 1,  Hunt  and  Hraham,  1978).  At  least  in  the  initial  stages  of  the  dis¬ 
tortion  (near  the  leading  edgei.the  horizontal  dimensions  of  the  eddies  just  above 
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the  viscous  layer  are  proportional  to  the  size  of  the  eddies  in  the  mainstream.  The 
author  (Lons,  1978)  has  discussed  the  importance  of  this  distortion  in  problems  of 
the  mixed  layer  in  stratified  fluids  in  laboratory  and  geophysical  situations.  In 
problems  of  shear  flow  there  is  a  large  production  of  eddy  energy  near  the  surface 
but  surely  the  large  eddies  of  the  outer  region  must  "feel"  the  wall  in  ways  similar 
to  those  of  shear- free  turbulence.  Townsend  (1976,  p.  161)  regards  the  turbulence 
as  being  composed  of  two  components,  one  called  a  "wall"  component  whose 
dimensions  scale  with  distance  from  the  wall,  and  the  other  an  "inactive"  com¬ 
ponent  with  horizontal  dimensions  of  the  order  of  the  outer  length,  and  vertical 
dimensions  of  the  order  of  distance  from  the  wall.  Perry  and  Abell  (1977)  call 
these  components  '  universal"  and  "non-universal",  respectively,  and  we  will 
adopt  this  latter  terminology.  The  two  components  of  motion  appear  clearly  in 
spectra  at  high  Reynolds  numbers  (Perry  and  Abell,  1975,  Bullock,  et  al. ,  1978). 

The  theory  in  this  paper  employs  independence  of  a  near  the  wall  and 
independence  of  v  far  from  the  wall  in  the  first  approximations  so  that  there  are 
still  the  sublayer  and  the  outer  layer.  It  appears,  however,  that  there  is  a  second, 
thicker  viscous  boundary  layer,  which  we  call  the  mesolayer,  whose  existence  must 
be  fundamental  for  the  problem  of  turbulent  shear  flow  in  general,  and  that  this  layer, 
intruding  between  inner  and  outer  layers,  prevents  the  overlap  assumed  in  classical 
theory!  The  mesolayer  appears  to  be  the  layer  in  which  the  horizontal  velocity  of 
the  large  eddies  of  the  non-universal  motion  reduces  to  zero  at  the  wall  as  in  the 
experiments  of  Uzkan  and  Reynolds  (1967).  Indeed,  we  find  that  its  thickness,  as 
determined  theoretically  in  Section  2,  is  identical  to  the  empirical  thickness  R‘ 
of  the  viscous  layer  found  by  Uzkan  and  Reynolds.  The  layer  is  similar  also  in 
eharacterand  behavior  to  the  laminar  boundary  layer  in  flow  above  an  oscillating  plane. 


1  The  failure  of  the  classical  approach  is  shown  in  detail  in  Section  4. 
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2.  Discussion  of  the  Mesolayer. 

In  Section  1,  we  presented  some  general  ideas  about  turbulent  shearing  flow. 
To  be  definite,  we  now  confine  attention  to  a  smooth  pipe  although  the  ideas  have  very 
general  applicability.  We  may  write  down  the  two  non-trivial,  exact,  mean  equations 
of  motion  (Laufer,  l954i. 


du 

dz 


uw  =  1  -  - 

R. 


WW  )  = 


ww  -  vv 

Rt(1-C) 


(5) 

(6) 


where  z  is  distance  from  the  wall,  so  that  z  R^  at  the  center  of  the  pipe,  C  =  z/R^  , 
p  is  the  mean  pressure  and  the  velocities  are  u,  v,  w  in  the  streamwise,  transverse 
and  radial  (inward)  directions.  As  in  the  rest  of  the  paper,  we  scale,  generally,  on 
inner  variables  v  and  u_ .  The  Reynolds  stress  -uw  is  zero  at  the  wall  and  its  first 
derivative  with  respect  to  z  is  zero  at  the  wall.  Moreover,  -uw  is  zero  at  the  axis 
of  the  pipe  so  that,  as  measurements  also  show,  it  is  a  maximum  at  some  distance 
from  the  wall,  say  z  -  m(R  ).  Therefore,  u22  =  -l/R„  both  at  z  =  0  and  z  =  m. 

When  R_  is  large,  the  effect  of  viscosity  is  confined  to  a  thin  layer  and  the  Reynolds 
stress  and  the  total  stress  (1-z/Rt)  are  nearly  equal,  except  at  very  small  z.  Con¬ 
sequently,  the  maximum  of  -uw  must  become  closer  and  closer  to  1  as  Rt  increases. 
Thus,  near  z  =  m,  we  have1  u2  —m/R^  . 

Let  us  now  make  the  tentative  assumption  that  u2  ~  m_1  in  the  region  of  the 
maximum  of  -uw.  The  two  order-of-magnitude  estimates,  u  ,  ~  m" 1  and  u2  —  m/Rr 

1  We  say  A(RJ  ~  B(Rt)  if  A/B  tends  to  a  finite,  non-zero  constant  as  RT-  *. 

We  also  use  A  =  0(B)  to  indicate  that  A/B  tends  to  a  constant  (independent  of  Rt)  which 
may  be  zero.  A  point  z  whose  distance  from  the  wall  varies  from  experiment  to  experi¬ 
ment  (i.e. ,  is  some  function  of  Rt)  is  said  to  lie  always  in  the  mesolayer  if,  for  all 
large  R  ,  z  =  om(R^)  where  a  is  a  constant. 
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lead  to  m  ~  R‘  .  We  have  an  indication  in  this  result  of  a  new  boundary  layer  over 

which  certain  quantities  change  from  their  universal  wall  values  to  those  typical  of 

the  outer  region!  This  second  layer  is  much  thicker  than  the  viscous  sublayer  over 

which  the  mean  velocity  changes  rapidly.  It  is  quite  evident  in  observations,  for 

example,  by  Laufer  (1954,  his  Fig.  8).  If,  for  the  moment,  we  use  the  empirical 

fact  that  u.j  is  very  close  to  l/n  7?  where  h  is  von  Karman's  constant,  we  predict 
_  1  _i 

the  maxima  of  -uw  at  /.  a4R*  where  a4  =  (h)  2,  or  at  values  of  C  of  0.  01 6  and 
0.  047  for  Laufer's  two  experiments  in  very  close  agreement  with  his  Figure  8. 

In  contrast,  the  maximum  of  ou ,  for  example,  is  very  much  closer  to  the  wall  and 
is  just  outside  of  (fixed  with  respect  to)  the  viscous  sublayer  (Laufer's  Fig.  25)  as 
we  discuss  further  below. 

Perry  and  Abell  (1975)  found  a  layer  (in  inner  variables)  in  which  the  data 
curves  of  ou  tend  to  become  level,  indicating  in  their  view  that  the  classical  theory 
correctly  predicts  o„  const  in  a  region  1  «z  «  Rt>  From-  our  present  viewpoint 
the  inflection  point  of  the  ou -profile  found  in  all  measurements  should  lie  in  the 

mesolayer.  Indeed,  these  points  are  further  out  from  the  wall  at  higher  Reynolds 

.1 

numbers  in  rough  proportion  to  R^'  as  we  will  also  see  below. 

T 

As  we  have  discussed  in  Section  1,  the  mesolayer  appears  to  be 
similar  to  that  in  shear-free  turbulent  flow  above  a  surface,  as  in  the  oxperi- 


lThe  classical  approach  also  predicts  a  maximum  of  -uw  at  7.  -  R^  but,  in  ignoring 
the  physical  effects  of  the  mesolayer,  it  fails  in  other  fundamental  ways  as  we  see 
in  Section  4. 


merit  of  Uzkan  ami  Reynolds  (19071  or,  perhaps,  in  turbulent  flow  near  a  surface 
at  some  distance  from  an  oscillating  grid  (lloptingor  and  Toly,  197M  as  suggested 
b>  the  author  (Long,  L9781.  The  present  situation  is  not  a  direct  analogy  to  the 
shear- free  problem,  however,  because  there  is  already  one  viscous  layer  at 
the  wall  of  the  pipe  and  it  is  not  obvious  how  to  obtain  another  one  trom  a  single 
set  of  equations.  We  may  obtain  both  as  follows,  however:  we  imagine  an  initial- 
value  problem  of  a  very  long  pipe  immersed  in  a  Navier-Stokes  fluid  initially 
at  rest  with  a  suddenly  imposed  force  along  it  at  t  0  creating  a  stress  if  on 
the  walls  which  is  then  held  fixed  in  time.  The  pipe  will  move  and  eventually 
come  up  to  a  constant  speed  u0>,  (a,  v,  u  ).  We  assume  complete  determinism' 
of  the  tuibulcnt  motion  such  that  the  velocity  component  q,  relative  to  the  wall, 
for  example,  is  the  same  at  (\,« yd, ttl  ;  a,  v  ,u  1  in  any  two  experiments  if 
the  initial  conditions  are  identical  even  after  a  very  long  time  when  the  problem 
becomes  statistically  steady.  At  a  fixed  /,  as  the  radius  of  the  pipe  a  gi'es  to 
infinity,  we  assume  that  the  effect  of  a  becomes  weaker  and  weaker  so  that 
q,  (\i . y,i t ,  ;  a,  v,  u  1  approaches  a  function  ullt  (\, ,  y.,,  ,  t,,;  v,  ul.  The 

assumption  that  the  limiting  function  ui(  exists  means  that  we  may  express 


Ux  +  ur 

(71 

Tix  +  Up 

(SI 

where  we  have  gone  over  to  non-dimensional  notation.  iq  and  are  independent 
of  U  and,  as  we  have  constructed  the  problem,  u  and  u0  tend  to  zero  as  H) 

'The  reader  who  is  accustomed  to  thinking  of  turbulent  flow  as  "random" 
may  read  an  amplified  discussion  In  Chcrn  and  Long  (198th. 
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by  lotting  a  >■ .  1'lion,  if  wo  substitute  into  tin-  sot  ot  Navior  Stokes  initiations, 
whiolt  \\(>  donoto  as  N.ami  lot  U  •  ,  wo  obtain  a  sot  of  equations  \,  involving 

only  u* ,  v*.  iv,,  p*  and  inner  variablos,  v,  z.  t.  1'lto  form  is  tho  same  as 
for  flow  ovor  an  infinito  piano  in  tbo  mfmito  half  spaoo  above  w  ith  tho  constant 
(dimensional!  momontnm  transfor,  u  .  and  above  tho  sublayer  (o.  g.  ,  in  tho 
mosolayeri  independence  of  v  loads  to  u,  k-'  fit  ?t.\,  o„,  A*.  ov ,  A., 

\,.  oto.  whore  A,  A,  are  univoisal  oonstants.  Sinoo  N,  does  not  involvo 
H  .  it  is  valid  for  all  H  and,  therefore,  wo  mav  subtract  N,  from  N  to  obtain  a  now 
sot  ot  equations  valid  for  all  K  ,  namely  N  ..  Obviously  tho  non  nnivorsal  oonniouont 

of  motion  may  be  idontifiod  with  tho  largo  soalo  motions  assooiatod  with  tho 
ooiv.  If  wi'  integrate  tho  munition  of  continuity  ovor  a  "ho\"  in  tho  mosolayor 
of  vortical  dimensions  of  ovdov  m  and  hovi;ontal  dimensions  l  oorrosi'onding 
to  tho  peak  of  tlu'  energy  spectrum  as  a  function  of  stroamwiso  wave  number 
k,  1,*'  ,  w  o  got  mis  values  ou  ■  -  >’..:l  ,m  ’  .  According  to  our  basic 

concept  of  tho  mosolayor,  tho  horizontal  eddy  velocity  ot  tho  largo  eddies  is  tho 
same  order  in  tho  mosolayor  as  in  tho  core.  Ising  independence  of  v  ,  we  have 
,i  ov  .  ~  i  in  the  core  and  in  the  mosolayor  -o  that  o.  .  ml  in  tho  mosolayor. 
In  tho  shear  free  o\porimonts.  tho  horizontal  od.fv  dimensions  in  tho  viscous  layer 
are  tho  same  order  as  in  the  region  outside.  Wo  assume  then  that  l  ,  H 
so  that  o.  .  mil" '  .  I  ho  thickness  of  tho  mosolayor  may  now  bo  obtained  by 


integrating  the  streamwise  momentum  equation  in  N2  twice  over  the  "box'  and 
equating  the  largest  viscous  term,  of  order  o^R*  ,  to  the  inertial  term1  of 
order  o?-,R?  ma.  We  find  a  thickness  m  R^  .  Thus,  we  obtain  the  same 
expression  R^  for  the  boundary- layer  thickness  and  we  have  identified  the 
mesolayer  with  the  "sloshing"  of  the  large-scale  eddies  along  the  surface. 

The  analysis  in  this  section  depends  on  the  two  assumptions  for  the 
mesolayer  that  uz  ~  m-1  and  that  the  axial  length  scale  is  of  order  Rt  for  the 
second  component  of  motion.  The  first  will  be  discussed  further  in  Section  3. 

The  second,  as  mentioned  earlier,  is  in  accordance  with  observations  by 
Tritton  (1967),  Mitchell  and  Hanratty  (1966),  Bullock,  et  al. ,  (1978)  and 
others.  These  measurements  are  especially  noteworthy  because  they  show 
that  over  a  considerable  range  of  Reynolds  numbers  the  horizontal  length  is 
not  only  large  near  the  wall  but  scales  best  on  the  outer  length.  The  observations 
of  Mitchell  and  Hanratty  were  made  remarkably  deep  within  the  viscous  layer 
indicating  that  the  large  eddies  are  exerting  an  influence  even  for  z  <  1! 


1  We  use  a  coordinate  system  moving  at  the  mean  speed  of  the  fluid  at 
the  center  of  the  box  and  assume  u ,  =  0(z_1 )  in  the  layer.  This  assumption  is 
verified  a  posteriori.  The  argument  involving  the  integrals  over  the  "box"  is 
given  in  more  detail  in  Chern  and  Long  (1980). 
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•b  The  Now  Theory  for  Mean  Quantities. 

Let  us  now  construct  a  new  theory  incorporating  the  mosolaycr  as  a 
fundamental  element.  Quito  generally,  we  may  write,  for  example, 

u  Tt,  (z)  fli,(z  ,  R^)  (>)) 

where  z  is  a  new  similarity  variable  z  z/m  ,  where  m  '  R^  is  the  thickness  of  the 
mesolayerl  If  z  constant,  the  point  z  is  in  the  mosolaycr.  We  assume,  as  a  funda¬ 
mental  property  of  the  layer, that  all  mean  quantities  are  of  the  same  order  at  all 
point s  of  the  mesolayer,  e.  g.  ,  u -Tl,  u,(z,  R^).  Then,  for  some  function  lvfll,) 

T 

so  that 


Ti  u,  (z)  +h0  (R^u*o(z)  +  u*i  (z,  R)  (10) 

"  hore  u4l  in  the  mesolayer  is  small  compared  to  lv(Rju,0(z).  The  same  argument 
tor  utl  shows  that  it  must  also  be  of  the  form  hi  (Rr)u+1(z).  Thus  we  may  write 


u  u,(z)+ho(R1Wi,0(j)+h1(RT)Il+1(3?)+...  (in 

T  T4(z)+h50(RT)l\0(7:)+h61(RT)T+1(zH...  (12) 

where  T  -uw  is  the  Reynolds  stress.  The  assumption  that  the  first  component 
of  motion  is  the  only  surviving  component  as  R  -  "  at  fixed  z  yields  two  equations 
in  place  of  (f») ,  namely 

»:  +  T,  1  (13) 

_  1  T 

li0R/u;o-ih1R/«;i  +  ...  +hr,0T*tffh„jT,1  +  ...  -  C  (14) 

where  primes  denote  derivatives.  Since  T  and  its  first  derivative  are  zero  at 
z  0,  we  have 

'  In  Section  I,  we  derive  m  -  R‘  by  arguments  independent  of  those 
of  Seel  ion  2. 
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bc  itu';o(0)--i,  u*i’(0)  0,  u;-?(0)  o . xs;x(0>  o,u;-<o>  o... 

Evaluating  (11)  in  t ho  mesolaver,  we  find 


»Vo  R,‘  ,  hx  h..lf  lyJl  *  h.,?, 


Solutions  for  Moan  Velocities,  Reynolds  Stress,  Dissipation  Function  and 
Energy  l’lux  -  Divergence . 

Lot  us  now  match  the  first  approximation  for  the  form  for  the  velocity  - 
defect  in  the  region  near  the  wall  with  its  form  in  the  outer  region  assuming 
independence  of  v  in  the  outer  region.  We  have ,  to  the  present  order, 

Uo  -  u  cp  (C  )  (Mo)  (15) 

u  u*(z)  +  Ti+0(z  )  (M, )  (ltd 

where  we  denote  here  and  below  regions  (0  <  z  «  R  )  by  M, ,  (R“  «  z<  ID  by 
(R‘  <•  /.  ■  llT)  by  M,  0.  Using  ^  u0(R  ),  we  have 


u0(Rt)- u,  (/.)-Ti+0(z)  c r(C  )  (Ml0) 

Differentiating  (17)  with  respect  io  R  .holding  '  fixed,  we  get 


(17) 


R 


du,,  ,  du.  I  *  tfu*o  ...  v 


dR  dz 


(IS) 


Differentiating  (18)  with  respect  to  z,  holding  R  fixed,  we  get 


>r  —  ■/  ^  i  ;  0  (ft!  \ 


(it» 


dz  dz  2  dz  dz. 

The  terms  in  (19)  must  each  equal  a  constant  so  that  the  solutions  for  u,  and  u0 
are ,  to  the  present  order, 

u,  C0  in  z  +  C,  (z  •»!)  (20) 


Uo  A1;.+  Cfft\+(C0+lC00)tn  Rr 


(21) 


lo  and 


(z  ■»!) 
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Also 

u.o  C00fn  z  +  C,  (MJ0)  {2.1) 

Uo  -  u  Axa-!C0+CooUn  C  (M,0)  (23) 

where  C0,  ,  Air,  C\,  C00  are  universal  constants  and  we  reject  a  solution  of 

the  form  z~l  in  z  for  diu/dz  because  of  the  principle  of  independence  of  v. 

The  expression  in  (20)  for  the  mean  velocity  is  the  same  as  in  classical 
theory  although  they  have  different  regions  in  which  they  are  supposed  to  apply. 
Obviously,  C0  is  close  to  >~x  where  h  is  von  Ka'rman's  constant.  The  Coe,- 
terms  in  (21)  and  (2:?)  are  absent  in  classical  theory,  but  we  have  no  reason  for 
assuming  this  in  the  present  analysis,  even  though  experiment  indicates  that  C0o 
is  small  compared  to  (J,.  The  fact  that  the  profile  of  u  is  logarithmic  both  below 
and  above  the  center  of  the  mesolayer  with  the  contribution  of  the  second  component 
small  means  that  u  (as  opposed  to  <tu  ,  <rv ,  as  we  discuss  below)  is  a  poor  indicator 
of  the  mesolayer  and  explains  why  the  classical  theory  has  been  so  enduring. 

If  we  use  (20)  and  (22)  in  (Id)  and  (14),  we  find 


T* 

1-CcZ-1 

(a  »l) 

(24) 

T*o 

-z-CcrZ-1 

(MI01 

(25) 

Notice  that  the  equation  T/+  R"1  T;0  0,  or  Cora  + T;0(z)  0,  defining  the  position  of 

the  maximum  of  the  Reynolds  stress,yields  z  a,  ,  showing  that  the  peak  of  the 
Reynolds  stress  is  in  the  mesolayer  as  suggested  above  and  by  the  data. 

A  physical  interpretation  of  our  results  is  that  dTgVdz  represents  the 
shear  in  a  region  z  ^  1  above  an  infinite  flat  plate  and  since  a  is  absent  and 
viscosity  should  become  negligible,  zdfi+/dz  must  be  a  universal  constant  C0  on 
dimensional  grounds  in  agreement  with  (20).  When  R  is  finite,  zdu/dz  C0  should 
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hold  in  some  region  near  the  wall.  On  the  other  hand  in  M10(R*  <sc  z  «  R)  both  v 
and  a  are  unimportant  and  we  again  find  zu,  =  C0+C00  is  a  constant  in  agreement 
with  (23). 

Let  us  now  assume  strict  independence  of  v  in  the  outer  region  to  solve  for 


T  in  Mio-  We  have 

-i 

T++R/T*0  =Xo(C> 

(Mj  0) 

(26) 

Matching  yields 

T+  =  *1— *  jz'1 

(z  »1) 

(27) 

T  *0  =  +  Hg  Z  1 

(Ml0) 

(28) 

Since  x  =  1~C  to  this  order,  nx  =  1,  k3  =  -1  and  by  comparison  with  (24)  and  (25), 

-  "Coo  =  C0.  According  to  (23)  we  lose  the  logarithmic  behavior  for  the  velocity 
defect  in  Ml0.This  seems  highly  unlikely  but,  in  any  case,  the  deduction,  C0+C00  =  0, 
cannot  be  made  once  we  recognize  that  independence  of  v  cannot  be  strictly  correct. 
Thus  we  are  required  to  add  a  higher  approximation  on  the  right-hand  side  of  (26) 
of  the  form 

T++R^T%0  =  l-C+HT/xifC)  (M,0)  (29) 

If  we  now  match,  we  find  that  the  coefficients  of  z-1  and  zT1  in  (27)  and  (28)  need 
not  be  the  same  and  we  get  the  behavior  of  (24)  and  (25)  without  the  requirement 
that  Cgt  Cqq  —  0. 

We  may  make  two  interesting  interpretations  of  the  mesolayer  which  help 
to  overcome  the  difficulty  in  understanding  how  molecular  friction  can  be  important 
in  a  region  far  above  the  sublayer.  We  may  define  the  sublayer  as  the  region  in  which 
the  Reynolds  stresses  are  of  the  same  order  as  the  viscous  stresses.  In  the  meso¬ 


layer  the  viscous  stress  becomes  negligible,  as  we  would  expect,  i.  e. , 
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Ti,/-(u’w')  -R.  .  However,  -u'w’  is  nearly  constant  in  the  mesolayer  (indeed  the 

classical  logarithmic  layer  is  often  referred  to  as  the  "constant  stress”  layer)  so 
that  its  derivative  is  small:  (-u'w'l  ~  z  +  HT 1  Tl0  or  (  u'w').  ~R_1.  This  is 
exactly  of  the  order  of  u. ,  revealing  that  in  the  mesolayer  the  forces  duo  to  the  Reynolds 
stress  are  ot  tin*  same  order  as  the  forces  due  to  the  viscous  stresses  and  that  both 
are  very  small.  Another  interpretation  of  the  viscous  effects  in  the  mesolayer  is  to 
conjecture  that  the  "bursts"  originate  in  the  region  very  near  the  wall  and  then  grow 
by  molecular  diffusion,  with  dimensions  of  order  (vt)*.  In  a  time  period  of  the 
order  of  the  large  core  eddies  this  becomes  (H\/u  )*■,  i.  e. ,  of  the  thickness  of  the 
mesolayer.  Apparently  the  bursts  transport  most  of  the  momentum  and  the  small 
change  with  height  of  the  momentum  flux  is  controlled  by  the  viscous  growth  of  these 
bursts. 


Further  discussion  requires  use  of  the  energy  equation  (Uniter,  10M): 

(30) 


where  f  is  the  dissipation  function  non-dituensionalized  by  inner  variables. 


dll’,  chi’,  ,  l 

•V.  A- .  '*-%/* 


,  dw'  rtV 


Ox.;  Ox,  ,  T  R  (l  -z.R  ~  1 )  V’  ay  +  RJl-zR*') 
and  F  is  the  non-dimensional  energy  flux 


(M,)  (31) 


where  P’  is  the  turbulence  pressure  function,  q’  is  turbulence  speed,  u\  v’,  w’ 
or  u',  are  turbulence  velocities  and  all  quantities  are  scaled  on  u^  and  v  .  Writing 
f.  +  b*o  »  in  (Mt)  and  using  Kq.  (d4)  and 


dh,  dF+ 
♦  dz  *  dz 


(3.!) 


C0  7.“  1  + 


‘IF, 

dz 


(z  >'  1)  (03) 


w  e  get 
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Suppose  (IF.  (1/  —  /“ Then  F,  behaves  logarithmically  ami  therefore,  because  all 
correlation  coefficients  for  the  first  component  of  motion  must,  to  first  order,  be 
universal  constants  for  /.  '»'»  I,  velocities  in  the  mesolayer  arc  of  larger  order  than 
one,  contrary  to  our  fundamental  assumption  in  deriving  m  -  ltf  from  the  N .  set  ol 
Navier-Stokes  equations.  Also,  of  course,  F,  ~  in  /  violates  the  principle  of  inde- 
pendence  of  v  .  We  reject  this  logarithmic  behavior  and  require  that  F*  be  smaller. 
Then  it  follows  from  (331  that 

€»  C0  1 ,  F. 

where  IT,  is  a  new  universal  constant, 
get,  to  first  order, 

H-,0R/  (M10>  (35) 

This  shows  thiit  lv0  R*  and  l\v  1.  Applying  independence  of  v  to  F  in  the  outer 
region  anil  matching,  we  find  that  F „ 0  is  constant  in  Mjo*  Equating  coefficients  of 

l 

R,*  in  (351,  we  find 

e„o  Coo  z  (M,0> 

This  completes  the  search  for  first  approximations  for  e,  ut,  T,  and  F,  but  it  is  now 

clear  from  our  presont  results  that  a  continuing  process  is  possible  in  which 

1  ,  -i 

mean  quantities  have  the  form  e  r,  +  R  'C  g»i+...,  F  F.  +  F<0  +  R  * l'*i  +  •  •  • 
etc.  where  e+1 ,  Ftt  are  functions  of  ?,  and  t luit  we  may  go  on  to  higher  approximations. 

We  have  seen  so  far  that  T  is  of  the  form  T  \0 ( 7  1+  RT1  \X(C)  in  the  outer 
region  (Mo '.and  tho  energy  equation  indicates  that  c  !*<('  ' K”,  +  ('  'RTt  , 

that  S  /.Tl,  and  F  are  of  similar  forms  and,  in 


(7.  ^  1) 

If  we  subtract  (321  from  (301,  we 
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general,  that  all  four  quantities  are  expandable  in  infinite  series 

approximations  are  obtained  from  the  matching, 

in  R” 1  • 

The  next 

l 

+  S+0+ RT  S+x  + . . .  =  ctq+ cpj  RT1  + . . . 

(M10) 

(37) 

T+  +  R~3T+0+  RT1T*1  +  . . .  =1-C+Xi  «?+*«?  +  ... 

( M  j  o ) 

(38) 

1 

F.+  F*o+Rt2F  +  1+.  ..  =  Aq+  Aj  rt1  +... 

(M10) 

(39) 

e+  +  R~*c0  +  R-T1e+X  +  ...  =  L0RT1  +  L1K;3+... 

(Ml0) 

(40) 

Retaining  the  first  three  terms  on  the  left  of  (37)  and  two  terms  on  the  right  and 
cross-differentiating  to  eliminate  first  S*  and  then  S+0>  we  find  an  equation  of  the 
form 

1 

*rt8(s*i+s;xS)  =h0(c)  +  r;1h1(o  (m10>  (4d 

so  that 

S*i+zS;i  =  2C10z+K2z_1,  S*1=C10z+C11z-1+K2z“1£nz  (M10) 

where  C^,  Cu,  K?  are  constants.  Eq.  (37)  is  now  of  the  form 

S*  +  S+0+K?R^£-1  Aiz  =H0o(C>  +  (Mt0) 

Eliminating  S*  and  solving  for  S+0,  we  get 

S*o  =  K4  in  z  +  C01z"3  +  Co0  (M10) 

where  Kj,  C0i,  C0o  are  constants.  Then 

S+  =C0  +  C!Z-1,  K4  =  K?  =  0  (z  1) 

where  C0,  Cx  are  constants.  Continuing  in  this  way  for  all  the  quantities  in  (37)  and 
(39),  we  obtain 


S+  =  Cq  +  Ci  z- 1  +  Cs  z-3  + . . . 

(z  >>  1) 

(42) 

S+0  -  C0o+  Cq \7.  +  C^ ^z  + .  •  • 

{M|  r, 

(43) 

S*i  =  C10z  +  C11z"1  +  C1^z_3+... 

(Mio  1 

(44) 
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F+  =  D0+Diz”1+D2z~ii+... 

(z  »  1) 

(45) 

F+0  -Doo+D0iZ  +D02z  +... 

(Mj  o  ) 

(46) 

F*i  =  Di0z  +Dnr1  +  Di2z"3  +  . . . 

(M  i  o  ) 

(47^ 

If  we  multiply  Eq.  (40)  by  z,  we  obtain  by  an  identical  argument 

<?+  =  E0z-1  +  Exz  °+  E2z-3+  . . . 

(z  >>1) 

(48) 

e*0  -  E00z-1  +  E0iZ-3+... 

(M10) 

(49) 

€+i  =  Eio+  EnZ  F . . . 

(M10) 

(50) 

Similar  arguments  applied  to  (38)  yield 

T+  =  B0+BiZ  1+B2z~*i  +  ... 

(z  »1) 

(51) 

T+0  =  B00z  +BoiZ-i  +  B02Z-iJ+... 

(M,o) 

(52) 

T+i  =  Bio  +  Bnz  +... 

(M,0) 

(53) 

•  •  • 

where  (5),  (32)  and  the  difference  between  (30)  and  (32)  yield  a  number  of  relations 
among  the  constants  including 

E0  =  C0,  ~  -Cq  —  Dj,  ,  E00  —  Cqo  ,  Eq  i  =  Coi  — 2Dqi 

B0  =  i,  Bx  =  -C0,  B2  =  -Cj,  Bqo  -  -I,  B0i  -  -Coo,  Bq2  -  -Coi 
Bio  =  -Cio,  Bn  = -Cn»  Bg0  = -C20»  B2i  = -Cai,  ... 

In  outer  variables  we  get,  to  first  order. 


S  =  Cq+  Cqo+  Cio  £4"  Cgo  C  +••• 

(Mt0) 

(54) 

F  =  Dq  o  4"  Bq  +  Dx  o  C4"  Dgo  C  + . .  • 

(M,o) 

(55) 

eR  =  (C0o+C0)T1  +  Ei0+E?oC  + 

•  •  •  (Mi  o) 

(56) 

The  mean  velocity  field  is  given  by 

u0  =  Aig+C++C^C0+JCoo)AnRT+KiiR;4+KisR;1  +  ...  (57) 


_ , 
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Uo-Ti  Ai3-(C0+C0oWn  '-Cl0  '-*C30r-  +  -  • .  (M,0)  (58) 
u„  C. +  t\.  «n  z-Cjz"  ’  -jCrz  -*  C,z_  -. . .  (/  >'  1)  (59) 

An  attempt  was  made  to  eompute  the  eonstants  from  publisheil  data.  The 
solution  in  (57)  for  Uo  was  fitted  to  the  data  of  Perry  and  Abell  (1077)  to  obtain 
Alr  +  (\  +  t\  -4.1 55.  (C\  +  }C0;>)  2. -495,  K1X  2.  ;155.  The  comparison  is  shown 

in  Fig.  1.  Next  the  velocity -defect  law  in  (58)  was  fitted  to  the  data  of  Perry  and 
Abell  as  seen  in  Fig.  2.  This  yielded  Alr  1.50,  C\,  2.0:?,  (C,+  t\0)  2.  57 

so  tint  C\  2.  t>2,  Co0  -  0.25.  From  the  data  of  Laufer  (1954,  his  1  ig.  20),  we 
see  that  the  prediction  that  dF  dr  is  constant  agrees  well  with  measurements  and 
we  estimated  D^o  2.  We  now  have  the  constants  in  the  first  two  terms  in  (5o): 

eR,  1.66+  '' 1  (M10>  (oO) 

This  may  be  compared  to  the  estimations  of  Laufer  for  eR  .  For  the  higher 
Reynolds  number  and  at  smaller  '  the  agreement  is  excellent,  for  example 
25.  0  vs.  25.  50  at  '  0.  I,  15.  5  vs.  15.  51  at  C.  0.  2  and  9.  0  vs.  9.  5o  at  2  0.  5. 

The  measurements  of  eR  for  Laufer ’s  low  Reynolds  number  case  are  about 
15-20'  lower  than  these  three  predictions.  The  inner  solution  for  the  mean 
velocity  was  fitted  to  the  data  of  Laufer  to  yield  Ci  25.58,  C.  1.57.  The  com¬ 
parison  is  shown  in  Fig.  5.  We  now  estimate  D*  from  Laufer  (his  Fig.  19)  bv 
estimating  the  rather  small  difference  between  the  pressure  and  the  kinetic -energy 
diffusion  terms.  \’ery  roughly  we  get  Dj,  6.  We  no\  have  the  constants  in 
the  first  two  terms  of  (48): 

e,  2.  62r~ 1  +  12.  7z“"  (■/  T’"'  1?  (61) 

This  agrees  rather  well  with  the  data  for  e  in  Laufer’s  Fig.  19  for  z  '-20.  We 


» 
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eonehido  our  comparison  with  T.  in  (51).  The  result 

T*  1  -2.  62  z  1  -25.  58  z7~  (z  l)  (62) 

agrees  well  with  Laufer's  Fig.  26  for  z  >  10. 


Solutions  for  root -me:ui -spun  re  velocities  and  mean  pressure. 

We  have  found  estimates  to  first  order  of  rms  velocities  and  turbulence 
length  scales  in  the  mesolayer:  L*  ^  Ly  ~  IT ,  L,  ~R*,  ou2  ~avr  ~l,  r“‘. 

W'e  therefore  write 


tj(z)+  f+io+  hii  f*ii  +  hi?  f*i  j  +  •  •  • 

(M.) 

(63) 

=  fa(zl+  f*»+  h21f+?1+  h32f+sS  +  . . . 

(Mt ) 

(04) 

( ,(z)  +  Rt  f,3 q  +  hu  f+3 i  +  h3  a  f+x>  + . . . 

(M.) 

(05) 

where  h.  .  h, ,  ill.  )  and  f, ,,  f. (z  ).  W'e  also  need  Eq.  (6),  which  we  may 
w  rite,  to  first  order,  as  two  equations, 


RT‘(i-C)h40p;c>+ivtip;1+...+2R“sf3f;30+2f3’f^0+2ir1f+30f+'30 
+  2R?fJh31f,31  +2f3h31f;34  +  ... 


Kl  to  fa  -fef?+  2f3f+30R^+  ir1  f  J0+  2fs  h3l  f*31  + . . . 


(M, )  (66) 


-2fsf«o-f**o  ... 


PI  +  2f3f?’  -  0  (67) 

where  we  have  used 

T>  P*(z)+h40(Ri)P^0(z)+h,1(RT)piiI(£)+...  (Mt)  (08) 

Evaluating  (66)  in  the  mesolayer,  wo  get  k»0  r'-,  lx,!  R"1 . h3l  R"1, 

a?  R.  ,  ....  hxl  h31  R;J  h*?  hrP  R~‘  ...  .  In  the  outer  region  (1)^,.),  using 


<4 

kXo(0+  4iaiT»Flx(o+... 

(Mo) 

(<>•)) 

4 

14o(0+  (:i(Rt)F?x(0+... 

(Mll) 

(70) 

<r„ 

F..,o(C  )+  4i(Rt)F31(0+... 

(Mol 

(71) 

P 

F.,o(C)+4i(R  JF4X(D+... 

<M0) 

(72) 

wo  soo  that  K(j.  ((>)  boeonios 

( l  ~  ^  4«  >  ^-1  1  11"^  •  •  •  <-  ^  3 O  ^  30  *^30  ^  3  1  ^3  1  "t*  £  ^ 3  X  1*  3  0  ^  31  •  *  • 

(M ,)  (7;{) 

I4  ■<  Q+Zl:iol:Al  t;»i  t-  ^3 1  l4  3  i  t-  .  .  .  -  l4  33  ~*^14  3  3  t-j  ^  l4  3  ±  ^  I4  ^  ^  •  •  • 

so  that  tll  .e.,4  4i.  4  s  4  a  4a  4i  4i.  ... 

We  have  found  that  tho  energy  flux  in  the  region  z  ^  i  and  H  •"  is 

T 

*•’*  D0+  l)x/_1+  ....  Since  correlation  coefficients  in  this  region  must  be  of 

order  one  to  first  order,  we  get 

fflfj  const  +  const  z- 1  + . . .  (z  >~*  i)  (71) 

f.  On i  +  o ;i ? z~ 1  +  ...  (z  t)  (75) 


onst  z 

+  ...  (z  >^D 

(7-1) 

l+  ... 

<Z  *"'U 

(75) 

'+  ... 

(z  '  l ) 

(76) 

ic  outer 

forms  for  ou, 

arc  revealed.  We  match 


+  i7i..zT  1  +  o'i-, z  a  + . . .  +  f*l0(z)+ R  *‘f+ix(z)+. . . 

T 


<M,0>  (77) 


»'»o(0  +  k;xvix{C)  +  ... 

and  similar  expressions  for  ov,  <r„  and  p.  The  results  are 
"u  Oix+Pn+^pC  +Pi.-.ra+... 


"...  a'x x 6 °i s z  '  *ai a z  + 1 


<\  + bv?  r+  pp.-i  c  +• . . . 


(M  to  1  (VS) 
(Z'^l)  (70) 
(M,o)  (*0) 
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<b. 

Orl+fts~Z  ‘  f  Op.,  Z~  *■  F  .  .  . 

(*  >M1 

(811 

0.1 1 

+  S  j  +  St  3  C+  6  .,  T  C  +  •  •  • 

(8dl 

fry* 

“  1  1  — (i 

O  q  ^  t  (Vq  o  /.  ^  O  q  q  /*  »  •  •  • 

(z  >M) 

(8:11 

p 

04i  +  b4  ?r  +s4.,r3  +  --- 

(M|0) 

(811 

p. 

04  1  t  T  7.  +  .  .  . 

(z  *M) 

(85) 

The  constants  in  (78)-(S3)  wore  estimated  using  the  data  of  Perry  and  Abell  ( l i»7f>> 
for  the  outer  quantities  and  <v,and  I-auler  (19511  for  (rvt.  and  .  The  results  are 


oxi 

A  l  .o^;, 

A-1,,iXt  330 

(801 

Six 

0.13,  Sx? 

-1.534,  Bx;>  0.108 

(871 

f»r  i 

^”1.70,  tvs  ? 

-1.08,  -78 

(881 

Pa 

~  0  ,  8^ 

-1.48,  gL..-,  0.533 

(89) 

^  1 

1.07,  a,? 

-10.81,  52 

(901 

0  .  B,? 

-0.325,  S3.,  -0.0417 

(91) 

omparison  with  the  data  are 

shown  in  Figs.  4-9. 

As  indicated  by  Figs.  1,  5  and  10,  <ra  inereases  very  rapidly  to  a  peak  at 
/.  15  and  then  deereases  gradually  in  terms  of  inner  variables  and  sharply  in 

terms  of  outer  variables.  The  region  before  the  peak  nuiy  be 
regarded  roughly  as  the  region  in  which  u,  is  constant  so  that  if  we  represent  o„ 
as  of  order  zu,,  wo  get  the  strong  linear  increase  observed.  In  this  region  the 
production  term  in  the  energy  equation  which  produces  u-fluetuat ions  only  increases 
ou  and  has  a  minor  effect  on  <\.  and  cr*.  This  is  supported  by  the  data  (banter.  Fig.  1 
which  show  that  oj'  is  the  major  pail  of  the  total  turbulence  energy  in  this  region, 
lust  before  the  peak  of  iru,  however,  the  pressure-scrambling  terms  in  the  moment 
equations  (IFinze,  1975,  p.  .'Id 6)  begin  to  create  v-and  w-fluetuat ions,  ou  decreases 
and  nv,  it „  increase.  This  behavior  is  also  revealed  by  the  data  for  z  15.  The 


details  of  the  rms  velocity  profiles  in  and  near  the  mesolayer  are  of  great  importance 
but  we  need  first  a  derivation  of  the  behavior  of  mean  quantities  in  the  region  ■/.  ■  U* 
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Behavior  of  Mean  Quantities  in  the  Region  z  «  Rg  . 

For  computational  purposes  we  have  assumed  above  that  there  is  a  region 

(at  least  for  large  R^  )  in  which  U*,  au  +  ,  ...  are  good  approximations  to  u,  cru . 

and  this  is  a  good  assumption  because  in  experiments  there  are  regions  in  which 
there  is  good  independence  of  the  outer  length  (Reynolds  number  similarity). 
However,  as  we  increase  z  and  approach  the  mesolayer,  contributions  from  the 
mesolayer  arise  which  we  may  obtain  by  assuming  a  Taylor  series  expansion  of 
the  functions  of  z,  for  example,  f+10(z)  =  +  allz  z®+  ...  for  i  «  1.  We  get 

o\i  =  fx(z)+  umz  +  Q'lisz8  + . . .  +  R  ^(a^iZ*-’.-  ci2?za+ . . . )  +  ... 

b  (92> 

(0  <  z  «  R2) 


a„  =  fs(z)  +cySiiz  +<yaaz2+. . .  +RT2(®as1£  +  cy?2azi3  +  ...)+  ... 

h 

(0  <  z  «  R2  ) 


<TW  -  (z)  +  RT  (QJ3  pj  z^-  Cl/32  gZ  +•••)+  Rt  (<¥33  pZ  4*  O'^g^Z  +•••)  + 


(0  <  z  «  R‘  ) 


where  we  use  9u\,/9z  -  0  a'  z  =  9.  Also 


€  =  e+(z)  +  Rt2(Ex?0  +Eiaiz  +...)+  HT1  (Ej  30  + Ei31z  +...)+.. 

1 

(0  <  z  «  R2  ) 

—  T 


F  =  F*  (z)  *f  Dx  1 1  z  +  Dxx  ?z  +■•••+■  Ht  a  +#••)  +  •  •  < 

(0  <  z  «  R^  ) 


=  S+  (z)  -  Z3+C1U  z4*...  +  R'^(C1?4z4+C1p5zB+...)  +  . 


(0  <  z  «  R2  ) 


(97) 


T  l-z_1S+(z)-  R^  (VIU  z''+Cu  ./  F...)  -R~  (C'1?  z-'+C1?rz 


( ;  i  s  > 


whore  wo  hit vi'  used  (5)  and  required  1 1\:» t  T  V,  1'.,  0  ;it  z  0.  Also 

It  Tl*(z>  -  |za-*-1Cll4z4  f...  FR‘  :*r'-h  ..!F.. 


-  «rv  ir  i 

0  at  z  0. 

Also 

k  f'x P  j'  f.  . 

,  . )  F .  . 

(99) 


(0  •  7.  -  R“  1 


Tho  Behavior  of  Moan  Quantities  in  tho  Mosolayor. 

Wo  may  now  infor  tho  bohitvior  of  moan  quantities  in  tho  mosolayor  itsolt. 
According  to  theory  and  moasuromonts  in  shorn*  froo  turbulonoo  (Hunt  and 
Graham,  19781,  wo  expect  tho  mosolayor  to  tom!  to  cause  an  mow  asoflooal 
maxima)  in  tangential  volooitios  <ru  and  ov  sinoo  tho  wall,  in  offoot,  transfors 
energy  from  the  normal  component  o,  to  tlto  par:illol  components.  For  <>„ 
tho  inner  behavior  involves  an  additional  peak  at  fixed  z  near  tho  wall  as  wo  have 
already  soon.  From  (92)  wo  obtain 

,ru  0\i+(YxsZ  F  .  .  .  +Oxil'/  F  iViiqZ  F  .  .  .  F  R^'o'i  .•  X  7.  F  .  .  .  (M<()  (10th 

in  a  region  beyond  tho  inner  maximum.  Obviously,  at  very  largo  11^  ,  there  will 
ho  a  subregion  in  M, ,  in  which  o„  is  nearly  constant  for  a  long  distance  between 
points  Tx  and  TP  in  Pig.  10.  To  tho  loft  of  T?  tho  universal  law  of  tho  wall  is  obeyed. 
On  tho  other  hand,  at  very  largo  z,  but  not  near  tho  center  of  tho  pipe,  tho  behavior 
is  Reynolds  -number  independent  and  is  given  by 

o„  on  F  t\ i  F  (lOlt 

as  also  shown  in  Fig.  10  to  the  right  of  r4.  Hut  nearer  the  mosolayor.tho  behavior  is 


- - - 
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l 

<^U  +  +  +Bi;.zU^  +  ...  (Ml0)  (102) 

Since  we  have  assumed,  on  physical  grounds,  a  peak  of  »>u  in  1h<'  mesolayer  as  wc 
move  to  smaller  z  along  the  linear  curve  of  (101),  the  z  ‘-term  in  (102)  begins  to 
be  of  importance  and  ,  ■  0.  The  slope  becomes  less  and  less  negative  and 

finally  becomes  zero  (and  u„  becomes  a  maximum  at  Tj,  as  wo  portray  in 
Fig.  10.  Still  closer  to  the  wall, there  is  a  region  whom  none  of  the  forms  for 

°u  in  (100)  (102)  apply.  As  we  increase  z  in  the  region  just  to  the  left  of  T? 
the  a,  z  ....  terms  in  (100)  begin  to  be  important  and  the  curve  turns  up  as  shown 
to  meet  the  portion  to  the  left  of  T4. 

We  may  find  the  various  transition  points  in  Fig.  10.  Obviously,  in  the 
region  between  1?  and  14  and  not  close  to  cither,  the  only  important  term  from  the 
inner  solution  is  the  constant  term  and,  from  the  outer  region,  f  4  ^  0  ( z  ).  Thus 

(Tu  on  +  f+io(£' 

and  the  data  should  collapse  when  plotted  against  z  .  We  see,  in  fact,  a  clear 
tendency  for  this  behavior  in  tin*  data  of  Perry  and  Abell  (11)7 fi)  in  Fig.  11  which 
is  plotted  against  £n  z.  The  inflection  point  T,  is  at  fixed  z  which  we  see  from 
Fig.  11  is  approximately  :i.  0.  At  T;,  the  weak  z~ ’-behavior  of  the  inner  solution 
is  matched  by  the  weak  /'-behavior  of  the  mesolayer  contribution  in  F,q.  (100),  i.o. 

/_l  ~  '■  or  z~ft4T  o'  Ihi'  location  of  T?.  At  T4  the  z"a-term  matches  the  R~^  - 
term  in  (102)  and  this  occurs  at  z  -  R^  or  z  -  It*.  The  locations  of  TP,  T.,  and 
1'4  show  why  the  collapse  of  the  data  in  Fig.  II  is  confined  to  such  a  small  region 
at  the  typical  Reynolds  number  of  the  experiments. 

Not  all  of  the  details  of  this  description  of  ou  in  and  near  the  mesolayer 
can  be  verified  by  experiment  because  of  errors  of  measurement  and  too- 
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low  Reynolds  numbers  but  the  data  of  Perry  and  Abell  in  Fig.  11  is  in  remarkable 
agreement  with  the  general  behavior  outlined.  The  data  for  uv  from  both  Laufer 
and  Perry  and  Abell  is  less  reliable  but  the  theory  indicates  that  a  similar  be¬ 
havior  occurs  but  with  the  peak  in  the  sublayer  missing  and  with  cvx  in  (SI > 
negative  as  we  decided  in  (88).  For  ou,  is  tentatively  evaluated  as  ~0. 18 

in  Kq.  (87).  For  (xv,  the  data  does  not  permit  accurate  evaluation  of  and  we  have 
set  0  in  Kq.  (89). 

? 

The  outer  maximum  of  <ju  is  at  z  ~  II'  but  for  the  typical  Reynolds 

number  of  experiments  this  is  probably  indistinguishable  from  the  mesolayer 

region  /.  ~  R‘.  This  may  also  be  true  of  <rv  .  Very  roughly,  lender’s  data  indicate  a  peak 

of  o¥at  /,  a.. R‘  where  aP  9.0.  The  normal  component  <rw  has  a  maximum 

i 

at  /.  aa  R ‘  and  the  data  of  La u for  and  Perry  and  Abell  suggest  a-,  ~  7.  0.  This  is 
as  it  should  be  because  u,  is  little  influenced  by  the  second  component  of  motion  in 
the  inner  portions  of  the  mesolayer.  Because  of  this,  we  would  also  predict  that 
<rK  should  show  independence  of  R^  over  a  greater  distance  of  the  wall  region 
that  a,,  and  <xv  and  this  is  also  revealed  by  Laufer’s  data. 

The  analysis  of  u  in  the  transition  is  delicate  because,  as  observations 
show,  u  is  so  close  to  logarithmic  in  a  large  region.  We  may  arrive  at  some 
interesting  predictions,  however,  as  follows:  We  write  down  the  forms  for  S 
(which  is  the  slope  of  u  in  a  logarithmic  plot)  in  M,0and  M,,  as 

S  C0  +  C0o  +  C01£-a  +  ...  (M,0)  (109) 

S  C0-za+... 


(M„)  (l  01) 


We  also  need  the  outer  form  of  T: 


T  1  -i  U"-  -<C0+  Cpo  >  7.  ~x  R-s  -  Cc 


A  a.3  r«  ' 

■Q,  7.  R  +  ,  ,  . 


(Ml0>  (105) 


In  (1 0:>i  the  z  R.‘  -  term  probably  contributes  negative  curvature  so  that 
c  oi  -0.  A  schematic  plot  of  S  near  the  transition  region  is  shown  in  Fig.  12. 

The  upper  dashed  line  is  the  classical  solution  S  C0.  For  larger  z  and  very 
large  RT,  S  is  close  to  C0  +  C00but  at  smaller  z  the  C0i-term  causes  S  to 
increase  as  shown.  The  joining  requires  that  C00  is  negative,  perhaps,  C'p0  -0.  25 

as  calculated  roughly  from  the  data.  This  behavior  suggests  a  picture  of  u  vs. 
in  z  or  In  z  in  and  on  both  sides  of  the  mesolayer  as  in  the  dashed  curve  of 
Fig.  15.  The  data  of  Perry  and  Abell  (13751  is  suggestive  of  this  but  the  data 
points  are  so  crowded  that  one  can’t  be  sure.  Laufer’s  data  (Rotta,  1962)  also  has 
the  same  trends  but  again  the  deviation  from  logarithmic  behavior  is  weak.  The 
present  theory  does  not  apply  directly  to  turbulent  boundary  layers  at  zero  pressure 
gradient  but  a  similar  theory  in  preparation  (bong  and  Chen,  1950>  has  a  similar  behavior  for 
u  in  and  on  both  sides  of  the  mesolayer.  Feda  and  Minze  (1975)  measured  u  at  two 
Reynolds  numbers  and  both  am  suggestive  of  the  trends  in  Fig.  1.1.  Similar  measure¬ 
ments  of  Smith  and  Walker  (1959)  appear  conclusive.  Details  are  contained  in  Long 
and  Chen  (1980). 
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4.  Comparisons  of  the  Classical  and  Mesolayer  Theories. 

A  number  of  authors  n'onnekes,  1968,  Bush  and  Fendell,  1972,  Afzad  and 
Yajnik,  1973,  Afzal,  1976)  have  gone  on  to  higher  approximations  using  the  classical 
results  for  the  first  approximation.  If  one  assumes^Afzal,  1976) 

S  =  f0*(z)+  hj  (R^fJ (z)+  hgfR^fljfz)  +  ... 

T  =  go(z)+4(RT)gr(z)+£?(RT)g*(z)+  ... 
for  the  inner  expansion  and 

S  =  cp0(Q+mi(RT)cPi(C)+mg(RT)^?(C  )+  . . . 

T  =*0(0+  nl<RT)t1<Q+na<RT)t8(C  )+••• 

for  the  outer  expansion,  then,  substituting  into  (5)  we  see  that  h,  =  £,  =  m, 
or 

s  =  fo  +  ir1C+  R"2f/+  ... 

T  =gb*+  tc'gi**  irsg|+  ... 

and 

s  =  gb  +  K^1  %  +  FTJ.2  9s  +  . . . 

t  =  *o  +  k;1  *i  +  IT%  +  •  •  • 

Matching  yields 

S  =  C0t  Aojz”1 +A0sZ-,e+. . . +  RT1  (AioZ+An  + Alsz-1  + . . . ) 

+  R~2  (A30z“+AsxZ+As?+. .  • )  +  ... 

1The  forms  in  (106) -(109)  follow  from  the  universal  law  of  the  wall  and  Reynolds  number 
similarity  and  the  assumptions  that  in  regions  z  ~  1  and  C  ~  1  mean  quantities  have  the  same 

order  of  magnitude  at  each  point.  Afzal's  work  seems  to  be  the  most  recent  of  these  efforts 
and  he  claims  to  include  earlier  work  as  special  cases.  One  may  search  for  other  routes 
to  higher  approximations  in  hopes  of  rescuing  the  classical  first  approximation.  I  don't 

know  how  to  make  an  exhaustive  search,  however,  and  I  will  be  content  in  this  paper  to 
advance  the  new  theory  as  best  I  can  and, in  the  arguments  of  this  section,to  transfer  the 
burden  of  proof  to  anyone  still  favoring  the  earlier  theory. 


(106) 

(107) 

(108) 

(109) 
=  n,  =  R-1 

1  T 

(110) 
(111) 

(112) 

(113) 

(114) 
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We  may,  in  fact,  use  the  concept  of  the  mesolayer  as  a  layer  of  transition  from 
inner  to  outer  behavior  for  the  Reynolds  stress  T  to  derive  the 
thlokness  of  the  mesolayer  without  reference  to  the  "eddy -sloshing"  argument  of 
Section  2.  We  accept  as  sufficiently  general  the  following  forms  of  the  shear- 


function  S  and  Reynolds  stress  T,  in  which  the  similarity  variable  z  is  now  z  =  z/Rb 

T 


and  M, ,  M0  are  now  0  <  z  «  r  and  Rb  «  z  <  R  : 

—  T  T  —  T 

s=S+(z)  +  f  R"S‘S+1(z) 
i  =1 

(Mj ) 

(life) 

S  =  ctb(C  )+  r  R"a,cDt(0 
i=l 

(Mo) 

(117) 

T  =T+(z)  +  E  R'^T^fz  ) 
i  =1 

(M,) 

(118) 

T  =  Xo(C)+  T  R"blx,(C) 
i  =1 

(Mo) 

(119) 

where  aI  +  1  >s,,  a1+l  >  a, ,  ...  .  Using  a  =  (1-b)  1  ,  matching  yields 


=  A0+  E  A,z-a’  ,  T+  =B0+  T.  BlZ'b‘ 
i=l  i=l 

(120) 

=  A0-t-  E  A+10CaSl,  xo  =B0+S  B+I0Catl 
i=l  i=l 

(121) 

=  A+10  z  aSt  +  E  A  +  13za(S,"a]) 
j=l 

(122) 

=  ?  A+jlCa(SJ"ai)  +  A^-*1 
j=l 

(123) 

i  =B*10  zatl+  E  B*1Jfa(t*"bi) 
j=l 

(124) 

X,  =E  +  B,C"b! 


(125) 


Substitution  into  Eq.  (5)  yields 


bx  1.  I>tti  l+at,  ti  l-b,  sx  l-2b,  t,  s,+b 


Wo  obtain,  for  example,  for  larger  z 


-b- 


-l4-')  -  -  -ab 


1’  l-Coz-l  +  B?z  “  +  ...  +R  (-z+B„  nz 


+  11 


-b-s  . 


'(B«»  za(Sr'*b"D  +  ...l  +... 


(1J » 0 


(127i 


and  for  smaller  z  , 


T  L-C0z_1+Br 


+  ...  +  K 


l+b 


+. . . I  +. 


(1.2th 


Wo  find  in  either  case  that  the  existence  of  a  transition  region  z  -  1  requires 
i 


b  t  or  z  ~  R.  as  in  the  "eddy -sloshing"  argument  of  Section  2. 

We  may  now  conveniently  discuss  the  third  basic  assumption  of  classical 
theory  (Tennekes  and  Lumlcy,  1072,  p.  1-16  and  p.  2 bat  that  there  is  a  finite 
region  of  overlap  of  the  inner  and  outer  regions.  According  to  Tennekes  and 
Lumley  and,  indeed,  the  basic  concept  of  an  overlap  region,  v  should  first  be¬ 
come  important  at  large  11  at  some  z  ~  l  as  we  move  toward  the  wall,  and  H 
should  first  become  important  at  some  z  ~  II  as  we  move  outward.  As  portrayed 
in  Fig.  11,  the  overlap  region  at  high  It  is  the  double-hatched  region.  In  the 
present  theory  a  transition  region  replaces  the  region  of  overlap.  Here  both 
vand  II  are  important  but  H  decreases  in  importance  as  z  decivases  within  the 
region  and  v  decreases  in  importance  as  z  increases  in  the  region.  The  situation 
is  then  as  pictured  in  Fig.  1.1. 

Malkus  (197th  has  stated  that  an  overlap  ">*  not  required  by  either  logic  or 
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dynamics1'.  Nevertheless,  our  finding  that  viscosity  has  importance  in  a  layer 

much  great  in  thickness  than  the  sublayer  needs  explanation.  We  may  extend  an 

argument  first  advanced  to  me  by  my  student,  C.S."  Chern,  who,  together  with  the 

present  author,  extended  the  mesolayer  concept  to  the  problem  of  thermal  convection 

over  a  hot  surface  (Chern  and  Long,  1980).  In  this  case,  we  find  a  mesolayer 

thickness  6,  ~  II  !q  "  ,  where  K  is  the  coefficient  of  molecular  diffusion,  H  is 

the  depth  of  the  fluid  and  q  is  the  buoyancy  flux.  Chern  argued  that  a  thermal 

rising  from  the  hot  surface  would  move  vertically  but  also  be  carried  horizontally 

by  the  large  eddies  of  size  H.  Since  these  have  a  speed  of  order  (qH) 3  ,  their 
\  -if 

time  scale  is  T„  ll  q  "  and,  if  we  assume  that  the  thermals  expand  by  thermal 

A 

conduction,  they  attain  a  size  (KT,)2,  which  precisely  equals  the  mesolayer 
thickness  5..  After  that  they  presumably  lose  their  identity.  Thus  the  transport 
of  heat  by  these  thermals  which  must  be  of  primary  importance  in  the  layer  where 
they  exist,  i.e.  the  mesolayer,  is  influenced  directly  by  molecular  conduction 
because  they  increase  in  size  by  molecular  conduction,  even  when  they  are  far 
above  the  sublayer.  In  shear  turbulence,  observations  indicate  the  existence  of 
a  burst  phenomenon,  as  we  have  already  discussed,  and  the  bursts  may  play  a  role 
similar  to  that  of  the  thermals.  Then,  if  they  also  diffuse  by  molecular  processes, 

A 

their  size  after  a  time  of  the  order  of  the  outer  eddy  time,  H/ut  ,  is  (vfl/u^  )2  , 
which  is  again  precisely  the  mesolayer  thickness.  Such  may  be  the  basic 
importance  of  viscosity  in  the  mesolayer. 
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5.  Conjectures  on  Intermittoncy. 

Wo  have  conjectured  that  the  viscous  mosolayor  is  the  region  of  origin 
of  the  bursts  observed  in  pipe  and  boundary  layer  flows.  This  suggests  t  he 
existence  of  three  time-scales  in  this  layer.  One  is  the  time-scale  of  the  energy- 
containing  eddies.  Another  is  the  period  of  the  bursts.  This  may  be  identified, 
perhaps,  with  the  set  of  equations  Nr.  We  may  estimate 


chi'  u-  dii;.  ua 
.X  '  Tb  *“  Ul  i>x  '  R, 


yielding  T b  ~  U.  as  conjectured  by  a  number  of  experimenters  based  on  oh  so rvat ions 
(Rao,  et  al. .  1971,  Ueda  and  Ilin/.e,  1975,  HeidrieU,  et  al. ,  1977,  Chen  and 
lllaekwelder,  19781.  A  thin!  time  scale  is  the  duration  of  the  bursts,  Td .  The 
intermittoncy  of  the  turbulence  may  be  defined  as  v  T,  1\.  All  observationn  indi¬ 
cate  that  the  Reynolds  stress  is  provided  in  major  part  by  the  bursts  which  pri 
sumably  act  as  mixing  elements  so  that  we  may  write,  in  the  language  of  mixing- 
length  theory. 


mv  '  1  ~  (X  ly  A.  lb.  (1301 

where  %  1.,  u.  is  the  eddy  viscosity  and  we  have  written  lv  for  the  mixing  length 
and  u.  for  the  velocity  of  the  burst.  Kstimating  u.  -  1  and  writing  u.  ~  1  1C 


we  find 
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v~  IV  I, 


(13*  1 


If  wo  now  assume,  tiuito  naturally,  that  the  duration  ot  the  burst  is  the  time 

lor  travel  over  the  mixing  length,  we  got  Tj  ~  L. .  Hut 

v  Td  L,  J* 

T,  ~  Rt  L, 

so  that 


L,  -  R4  ,  v  -R/,  Td 


IV 


(13*5 


(133* 


The  expression  for  %  yields  \  0.  I  -  0.  2  for  typical  Reynolds  numbers  and 

this  is  in  rough  agreement  with  observations. 
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